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Abstract 

The quantum correlations between the beams generated by polariton pair scattering in a semi- 
conductor microcavity above the parametric oscillation threshold are computed analytically. The 
influence of various parameters like the cavity-exciton detuning, the intensity mismatch between 
the signal and idler beams and the amount of spurious noise is analyzed. We show that very strong 
quantum correlations between the signal and idler polaritons can be achieved. The quantum effects 
on the outgoing light fields are strongly reduced due to the large mismatch in the coupling of the 
signal and idler polaritons to the external photons. 

PACS numbers: 71.35.Gg, 71.36. +c, 42.70.Nq, 42.50.-p 
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INTRODUCTION 



Exciton polaritons are the normal modes of the strong light-matter coupling in semi- 
conductor microcavities [lj. These half-exciton, half-photon particles present large optical 
nonlinearities coming from the Coulomb interactions between the exciton components. Un- 
der resonant pumping, this leads to a parametric process where a pair of pump polaritons 
scatter into nondegenerate signal and idler modes while conserving energy and momentum. 
The scattering is particularly strong in microcavities because the unusual shape of the po- 
lariton dispersion makes it possible for the pump, signal and idler to be on resonance at the 
same time (see Fig. [TJ. Moreover, the relationship between the in-plane momentum of each 
polariton mode and the direction of the external photon to which it couples j2| enables to 
investigate the parametric scattering using measurements at different angles to access the 
various modes. 

The first demonstration of parametric processes in semiconductor microcavities was per- 
formed by Savvidis et al. using ultrafast pump-probe measurements. He observed para- 
metric amplification, where the scattering is stimulated by excitation of the signal mode with 
a weak probe field. Parametric oscillation, where there is no probe and a coherent population 
in the signal and idler modes appears spontaneously, has since been observed by Stevenson 
et al. Q] and Baumberg et al. in cw experiments. The lower polariton was pumped 
resonantly at the "magic" angle of about 16°. Above a threshold pump intensity, strong 
signal and idler beams were observed at about 0° and 35°, without any probe stimulation. 
The coherence of these beams was demonstrated by a significant spectral narrowing. 

The large optical nonlinearity of cavity polaritons makes them very attractive for quantum 
optics. Noise reduction on the reflected light field has been predicted (f| and achieved 
experimentally [7( for a resonant pumping of the lower polariton at 0°. The parametric 
fluorescence was recently shown to produce strongly correlated pairs of signal and idler 
polaritons, yielding a two- mode squeezed state The parametric oscillation regime is 
also very interesting in this respect It is well known that optical parametric oscillators 
(OPO) can be used to generate twin beams, the fluctuations of which are correlated at the 
quantum level. A noise reduction of 86% was obtained by substracting the intensities of the 
signal and idler beams produced by a LiNbOa OPO [lo| . 

The purpose of this paper is to investigate the possibility of generating twin beams using 
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a semiconductor microcavitv above the parametric oscillation threshold. The classical model 
developed by Whittaker [ll| is no longer sufficient to study the quantum noise properties 
of the system. Thus we adapt the quantum model by Ciuti et al, previously used in 



the context of parametric amplification [12] and parametric fluorescence |8|, [l3j, to the 



parametric oscillator confi gur ation. Furthermore we compute the field fluctuations using 
the input-output method llj, lal. We also include the excess noise associated with the 

nn 

excitonic relaxation, which had not been done in Refs. yj [13J . 



MODEL 

Hamiltonian 

Following Ciuti et al. |12j [13j we write the effective Hamiltonian for the coupled exciton- 
photon system. The spin degree of freedom is neglected. 

H — Hq + H exc _ exc + H sat (1) 
The first term is the linear Hamiltonian for excitons and cavity photons 



k k 

+ m R (4^k + &£a k ) (2) 

k 

with b k and a\. the creation operators respectively for excitons and photons of in-plane 
wave vector k, which satisfy boson commutation rules. E exc (k) and E cav (k) are the energy 
dispersions for exciton and cavity mode. The last term represents the linear coupling between 
exciton and cavity photon which causes the vacuum Rabi splitting 2MIr. The fermionic 
nature of electrons and holes causes a deviation of the excitons from bosonic behavior, which 
is accounted for through an effective exciton-exciton interaction and exciton saturation. The 
exciton-exciton interaction term writes 

H exc - exc = \j2 ^k+q^k'-q^k' (3) 

k,k',q 
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where V q ~ Vq = e ^ c for qa exc <C 1, a exc being the two-dimensional exciton Bohr 
radius, eo the dielectric constant of the quantum well and A the macroscopic quantization 
area. The saturation term in the light-exciton coupling is 

H S at ^ V sat (oi+q&ic'-q 6 * 6 * + ak+q&k'-q&M') ( 4 ) 

k,k',q 

where V sat = with 

n sat — 7/ (167rag XC ) being the exciton saturation density. We 
consider resonant or quasi-resonant excitation of the lower polariton branch by a quasi- 
monochromatic laser field of frequency u>l = El/H and wave vector ki,. If the pump intensity 
is not too high the resonances (i.e. the polariton states) are not modified. Then it is much 
more convenient to work directly in the polariton basis. It is possible to neglect nonlinear 
contributions related to the upper branch and consider only the lower polariton states. The 
polariton operators are obtained by a unitary transformation of the exciton and photon 
operators: 

f P ,\(-c k x k \ U\ 

\ <?k J \ x k C k J \b k J 

where X k et C k are positive real numbers called the Hopfield coefficients, given by 



. 2 s k + y/5 2 k + n 2 R 



X 2 = - ' v "ft ' (6) 

n 2 

C l = R N (7) 

2^[TWM k + ^i + W R 



X 2 and C\ can be interpreted respectively as the exciton and photon fraction of the lower 
polariton p k . In terms of the lower polariton operators, the Hamiltonian ([!} reads 

H = H P + H e Jl (8) 
Hp is the free evolution term for the lower polariton: 

H P = Y,E P {k)p{ P]i (9) 

k 

and H e Jp is an effective polariton-polariton interaction: 
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jjeff 

Upp 2 



\ ^SqPk+qPk'-qPkPk' (10) 



k,k',q 

where 



^Ck'.q = {VoX\k+q\X k > + 2V sat (11) 

X (C|k+q|^fc' + Cfc'^|k+q|) } -X"|k'-q|-Xfc 

In the following we neglect the contribution of the saturation term, so that V££ — 
VoX|k+q|^fc'X|k'-q|Xfc. We also neglect multiple diffusions i.e. interaction between modes 
other than the pump mode. This approximation is valid only slightly above the parametric 
oscillation threshold It comes to considering only the terms where the pump polariton 
operator pk L appears at least twice: 



H e Jl = - Vk L ,k L ,oP kL Pk L Pk L Pk L 

+ y y k L ,k L ,k L - k (pk L -kP k Pk L Pk L + h.c.^j 

k^k L 

+ 2 £ KIAAp^ ( 12 ) 

k^k L 

The first term is a Kerr-like term for the polaritons in the pump mode. The second term 
is a "fission" process, where two polaritons of wave vector kx are converted into a "signal" 
polariton of wave vector k and an "idler" polariton of wave vector 2kL — k. The last term 
corresponds to the interaction of the pump mode kL with all the other k states, which results 
in a blueshift proportional to |pk L | 2 - 



Energy conservation 

The energy conservation for the fission process {k L , k L } — > {k, 2k L — k} reads 

E P (k) + E P (2k L - k) = 2E P (k L ) (13) 

where Ep (q) is the energy of the polariton of wave vector q, renormalized by the inter- 
action with the pump polaritons 
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Ep(q) = Ep(g) + 2V q , kL , |(p kL )| 2 (14) 

Note that the factor of 2 disappears for q = kL. Equation (|13jl always has a trivial 
solution k = Icl- Non-trivial solutions exist provided the wave vector kL is above a critical 
value, or equivalently if the angle of incidence is above the so-called "magic angle" 8 C 0|. 
From now on we suppose that the microcavity is excited resonantly with an angle 8 C . Fig. |21 
is a plot of the quantity \E P (k) + E P (2k L — k) — 2E P (k L )| as a function of k = {k x , k y }, 
ki, being parallel to the x axis. 

This shows that energy conservation can be satisfied for a wide range of wave vectors 
{k, 2kL — k}. In recent experiments, parametric oscillation was observed in the normal 
direction k = jj, [l^. In this paper we consider only the parametric process {kL, kL} — ► 
{0, 2kL} assuming that the other ones remain below threshold. Then we can neglect the 
effect of modes other than 0, kL, 2kL- The evolution of these three modes is given by a 
closed set of equations. 



HEISENBERG-L ANGEVIN EQUATIONS 

In order to study the quantum fluctuations we have to write the Heisenberg-Langevin 
equations including the relaxation and fluctuation terms. The relaxation of the cavity mode 
comes from the interaction with the external electromagnetic field through the Hamiltonian 




The coupling constant is given by k = a/2t^ where 7 a fc is the cavity linewidth (HWHM). 
This leads to the following evolution equation for the cavity field: 

^(t) = - lak a^{t) + y/n^{t) (16) 

where A™(t) is the incoming external field. In this equation the normalization are not the 
same for the cavity field as for the external field: n ak (t) = /a^.(t)ak(t)^ is the mean number 
of cavity photons, while J^ 1 = (^A^\t)A™(t)\ is the mean number of incident photons per 
second. 



6 



Exciton relaxation is a much more complex problem. The density is assumed to be low 
enough to neglect the relaxation due to exciton-exciton interaction [l^j]. At low density 
and low enough temperature the main relaxation mechanism is the interaction with acoustic 
phonons. A given exciton mode is coupled to all the other exciton modes by and to all the 
phonon modes fulfilling the condition of energy and wavevector conservatio n UM- Relaxation 
in microcavities in the strong coupling regime has been studied in detail |l9j,|20|. However, 
the derivation of the corresponding fluctuation terms requires additional hypotheses, under 
which one can replace the exciton-phonon coupling Hamiltonian by a linear coupling to a 
single reservoir |9j. Then, in the same way as for the photon field, the fluctuation-dissipation 
part in the Langevin equation for the excitons writes 



-IbMt) + y/2 lhk B^{t) (17) 

where 7^ is the exciton linewidth (HWHM) and B^(t) the input excitonic field, which 
is a linear combination of the reservoir modes. 

Using these results we can write the Heisenberg-Langevin equations for the cavity and 
exciton modes of wave vectors 0, kL, and then for the three corresponding lower polariton 
modes. We define the slowly varying operators 

p (t) = Po (ty E ^ h 

p 2kL (t) = p 2 ^(t)e lEA2kL)t/h (18) 
which obey the following equations: 



° * (sV^k^opj^PkL - ijo) Po (19) 



dt h 

i 



^ ,k L ,k L P2k L Pk L e + r 



dp 2 k L i 



^2k L 

dt ~ h 

i 



dpk L 
dt 



(2^2k L ,k L ,opI L Pk L - «72fe L ) P2k L (20) 

* V n 1 J&Et/h , pin 

■^V^MlMlPoP^ 6 +i 2k L 

l - (A L + V^^opLjy^ - i-i k ^j p kL (21) 



2i 
T 



V kM AjoP2 klj e- iAE « n + P™ 



where for any given wave vector q, P™ = —C q ■ s j2'j aq A^ + X q ^J1^ hq B™ is the polariton 
input field (which is a linear combination of the cavity and exciton input fields ; only the 
driving laser field A 1 ^ has a nonzero mean value), 7 g = Cg7 ag + X^^ is the polariton 
linewitdth ; Ax = E p {k L ) — E L is the laser detuning ; AE = E p (2ki) + E p (0) — 2E L is the 
energy mismatch. 

These equations extend the model developed by Ciuti et al. above threshold. The 
full treatment of the field fluctuations is included as well as the equation of motion of 
the pumped mode accounting for the pump depletion. They are valid only slightly above 
hreshold, because far above threshold when we can no longer neglect multiple scattering 
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221. 



This set of equations is similar to the evolution equations of a non-degenerate optical 
parametric oscillator (OPO) (2^. The non-linearity is of type x^ while in most OPOs it is of 



3- 



However, 



type ■ OPOs based on four- wave mixing have already been demonstrated 
let us stress that here the parametric process involves the excitations of a semiconductor 
material (i.e. polaritons) instead of photons. In the following we evaluate the potential 
applications of this new type of OPO in quantum optics. The hybrid nature of polaritons 
makes the treatment of quantum fluctuations more complicated, since we have to consider 
additional sources of noise (i.e. the luminescence of excitons). 

MEAN FIELDS ABOVE THRESHOLD 

To start with we have to compute the stationary state of the system. This comes to the 
calculation done by Whittaker in Ref. We neglect the renormalization effects due to 

the interaction with the pump mode, which allows to get analytical expressions. Moreover, 
we suppose that the angle of incidence is adjusted in order to satisfy the resonance condition 
AE = and that the pump laser is perfectly resonant (Ax=0). Equations (|19 |) -(|21 |) now 
write 

^ = - 7o po - iE in A^ + Po (22) 



-72k L p 2kL - iE mtV \f^ + P™ L (23) 



dt 

^ = -7k L Pk L - 2iE in Aj p 2kh + Pg (24) 
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where E int = Vk L ,k L ,k L /^- Let us recall that among the polariton input fields, only the 
photon part of P£™ corresponding to the pump laser field has a nonzero mean value. The 
excitonic input fields B™ correspond to the luminescence of the exciton modes and are 
incoherent fields with zero mean value. The stationary state is given by 



- TfeiPkL - 2i£mtPk L PoP2k L = CWV 2 7aA< L ( 2 5) 

-7oP - ^mtP2k L Pk L = ( 26 ) 

-72fe L P2k L + iEintPoPil = (27) 

For a non trivial solution to exist, the determinant of the last two equations must be 
zero: 

£L|Pk L | 4 -7o72 fci =0 (28) 
which gives the pump polariton population threshold 



W =-*-Ei (29) 



E 



int 



and the pump intensity threshold 

it (7o72fcJ 1/2 



Tin 

J k L ,S - 



2 



Tin 



^laCl L E int 

The signal and idler polariton populations are easily derived: 



(30) 



i 2 7fcz, /72fci / -.x / Q1 x 
bol = Tr^J—iv-l) (31) 

^int V 70 



2-C/int V 72fc L 

where a = \Jl^Jl^~ s is the pump parameter. We finally get the intensities of the signal 
and idler output light fields 



jout 
J 



TOUt 

J 2k L 



2 7aC 2 |Po| 



iaik L Cl Fm 



= 2 laC 2kL |p 2kL | = 



Eint V 7o 
| 2 _ lalk L Cl kL 



{a -I) 



(33) 



7o 

72fc £ 
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Above threshold, all the polaritons created by the pump are transferred to the signal 
and idler modes, so that the number of pump polaritons is clamped to a fixed value. This 
phenomenon called pump depletion is well-known in triply resonant OPOs. The si gnal and 
idler intensities grow like yl^- These results are in agreement with those of Ref . . 

Finally we study the ratio of the signal and idler output intensities, which is an important 
parameter in view of the analysis of the correlations between these two beams. It is given 
by the simple equation 

rout /"(2 
-'o 72fc L O 

rout r<2 1° v 

J 2k L 70<- 2 fc L 

We consider a typical III-V microcavity sample containing one quantum well, with a Rabi 
splitting 2hQji= 2.8 meV. At zero cavity-exciton detuning, one finds £^=1.15 10 4 cm -1 . The 
photon fractions of the signal and idler modes are respectively Cq=0.5 and C\ kL ~ 0.053. 
Assuming that they have equal linewidths the signal beam power should be about ten times 
that of the idler beam. It is possible to reduce this ratio a bit by increasing the cavity- 
exciton detuning, as can be seen in Fig El However, the oscillation threshold goes up. In 
the following, all the results will be given at zero detuning. 



FLUCTUATIONS 



Linearized evolution equations 

For any operator 0{t) we define a fluctuation operator 50(t) = 0(t) — (0(t)). In order 
to compute the fluctuations, we use the " semiclassical" linear input-output method, which 
consists in studying the transformation of the incident fluctuations by the system |l5j . It has 
been shown to be equivalent to a full quantum treatment. We linearize equations (|22|l - (j24j) 
in the vicinity of the working point po computed in the previous section. We obtain the 
following set of equations: 



= -7k L 5pk L - 2iE mt (p p2k L <^k L (35) 

+Pk L P2k L ^o +P kL Po^2k L ) + SPg 
= -7o<*Po - iE int (2p2k L Pk L ^k L (36) 



10 



-l2kJP2^ ~ iE int (2poPk L ^Pk L (37) 
pIM) + SP 2^ 



dt 

tjPo 

We can now inject the mean values of the fields j?k L) Po et P2k L that we have computed 
in the preceding section (equations (|29|) . (}3Tf and (}32j) ). 

First we have to choose the phases of the fields (this choice has no influence on the physics 
of the problem). We set the phase of the pump field A™ to zero. Then p kL is a positive real 
number. The equations ()26|) and (}2Tj) impose the same relationship between the signal ipo 
and idler <^ 2 k L phases: 

7T 

<A) + <^2k L = -- (38) 

whereas the relative phase (po — <£>2k L is a free parameter. We set p to be a real positive 
number (again, this choice is of no consequence regarding the physics of the problem). Then 
p 2 k L is a pure imaginary number. With these choices of phase, the evolution equations write 



dSp klj 
dt 

d5p 
dt 

dSp 2 k L 
dt 



-Ik, 



(5p^ + (a-l)5pl h 



V27^7o (a-l)5p (39) 



-t^2 lkLl2kL (a - l)^ 2kL + 6P& 



-i^W^M^ + *Pt (40) 

-72fc L ^2k L - iV 2 lk L l2k L (O- - l)5pk L 

-iVWtolSpl + SP&i (41) 



Thanks to these three equation and their conjugate equations we can calculate the output 
fluctuations of the pump, signal and idler fields as a function of the input fluctuations. 

Amplitude fluctuations 



In this paper we are mostly interested in the amplitude correlations between signal and 
idler. We will see that in the simple case where we neglect the renormalization effects it is 
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enough to solve a system of three equations. We define the real and imaginary parts of the 
polariton, photon and exciton fields 



A, = ~i (5p q " Sp{) (42) 

in(out) r pin(out) , r pin(out)\ 

Ctq q q 

pin(out) _ _^ ^pin(out) _ fopin(out)]^ 

^A,in(out) ^j^in(out) _|_ ^^m(oiit)f 

pA,in(out) ^ ^^m(out) ^^m(out)f^ 

^B,in(out) ^ pin(out) _|_ ^-Qin{out)] 

pB,in(out) _ _^ ^pin(out) _ ^ gin(out)\^ 

The mean fields pk L an d Po are real positive numbers, therefore a corresponds to am- 
plitude fluctuations and (5 to phase fluctuations. The mean field p 2 k L is a pure imaginary 
number, therefore —f3 corresponds to amplitude fluctuations and a to phase fluctuations. 
The evolution equations for the amplitude fluctuations write 



= -7fc^"k L - y / 2-f kL -f (a- l)a 



+ V / 27 fei 72 fci (^-l)/9 2kL +< (43) 

dao i - 

— = -7 a + V 2 7fc L 7o \o ~ l)«k L 

-a/TO^/W + < (44) 

G^2k L 



-v/T^T^o + ^ L (45) 



We get a set of three linear differential equations. Taking the Fourier transform we obtain 
in matrix notation 



1 7 kL a-itt v /2 7^7o (o- - 1) -^7^72/^ (o" - 1) ^ / «k L (A) ^ 

\ & kL (n) 



V 



-V /2 7fe i 7o(^- 1) 7o 
V 2 7fcL7 2 fct (o- - 1) v^oT^I 



l2k L - ito 



X 



k L 

m 



(46) 
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The inversion of the 3x3 matrix provides the amplitude fluctuations of the fields pk L , 
Po et p2k L as a function of the input fluctuations. It is easy to deduce the amplitude 
fluctuations out of the output light fields thanks to the input-output relationship for the 
cavity mirror A™* = ^2 / -f aq a <i — A™ and the relationship between the photon and polariton 
fields a q = — C q pq. 

<' ou * = -C gV /^a q -< m (47) 

Input fluctuations 

In this paragraph we study the noise sources in our system. A™ is the coherent pump 
laser field A 1 q and both other input fields are equal to the vacuum field. Therefore, 
the amplitude fluctuations of these three fields are equal to the vacuum fluctuations. The 
treatment of excitonic fluctuation is more complex. The amplitude noise spectra (normalized 
to the vacuum noise) of the three excitonic fields B^, B l § et are gi ven by 

S£ m m = 1 + 2n q for q = 0, k L , 2k L (48) 

where n q is the mean number of excitations in the reservoir which depends on the temper- 
ature and on the pump intensity. Since the reservoir is populated through phonon-assisted 
relaxation from the pump mode it is a reasonable assumption to take the reservoir occupation 
as proportional to the mean number of excitons in the pump mode: 

n cl = P\b (l \ 2 = (3X 2 q \p cl \ 2 (49) 

where (3 is a dimensionless constant which characterizes the efficacy of the relaxation pro- 
cess. This simple model accounts for the excess noise of the reflected light at low excitation 

□ 

intensity in a satisfactory way |9J. 
Noise spectra 

In fluctuation measurements the measured quantity is the noise spectrum. The noise 
spectrum So{Q) of an operator O is defined as the Fourier transform of the autocorrelation 
function Co (t, t')\ 
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S (Q) = J Co(r)e inT dT (50) 

with 

Co (t, t') = (0(t)0(t')) - (0(t))(0(t f )) = (50(t)50(t')) (51) 

The noise spectrum is related to the Fourier transform 80 (Q) of the fluctuations 50 (t) 
by the Wiener-Kinchine theorem: 

(50(Q)50(Q')) = 27r5(Q + n')S (ty (52) 

In the same way the correlation spectrum Soo'ity of two operators 0,0' is defined as 
the Fourier transform of the correlation function: 

Co <y (f , if) = (0(t)0'(t')) - (0(t))(0>(t>)) (53) 
The correlation spectrum is also related to the Fourier components of the fluctuations: 

(50{Q)50'(Q')) = 2n5(Q + Q')S o'(^) (54) 
The relevant quantity is the normalized correlation spectrum 

Coo>(V) = ^£^B= (55) 

y/S (to)So.(to) 

One has always |C| < 1. A nonzero value of C o'{^) indicates some level of correlation 
between the two measurements. 



RESULTS 



Fluctuations of the intracavity polariton fields 

First, in order to shed some light on the above-mentioned analogy with an OPO, we 
assume that all three polariton modes have the same linewidths. This is the case if the 
cavity and exciton linewidths are equal {j a k = jbk) and do not depend on k. We set 

7 = lk L = 70 = l2k L =la = lb- 
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After some straightforward algebra we get the amplitude fluctuations of the polariton 
fields 



with 



- 7v /2 (a - 1) (2 7 - iO) < (56) 
+ 7v / 2(a-l)(2 7 -^)^ L ) 

a ° = DW) (W2(^-l) (27 - <n) < (57) 
+ ( 7 2 (3cr - 2) - fi 2 - i 7 fi (a + l))a™ 
+ 7 ( 7 (a - 2) + ifl) /^J 

& kL (n) = ^ (-7 VW^i) (27 - <n) < 

+ 7 ( 7 (a - 2) + in) a™ (58) 
+(7 2 (3(T _ 2) _ q2 _ i^a + i))^) 



D (n) = 7 [8 7 2 ((7 - 1) - fi 2 (tr + 2)] + iQ [ 7 2 (4 - 6a) + fi 2 ] (59) 
Twin polaritons 

Let us now calculate the fluctuations of the difference of the signal and idler amplitudes. 
Let r be the normalized quantity 

r = -^=(a + (3 2k J (60) 

We find 



r(n) = (4 7 2 (a - 1) - Q 2 - in-fa) r m (61) 
avecr- = -L( a »» + /^J 

It is important to notice that r does not depend on the pump fluctuations, which cancel 
out when we make the difference. This property is at the origin of twin beams generation 
in OPOs. We get perfect noise suppression for f2 = and a — > 1. 
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In a degenerate or quasi-degenerate OPO the symmetry between signal and idler is con- 
served outside the cavity, because the two fields have the same frequency and are coupled 
in the same way to the external field through the losses of the cavity mirrors. In such 
systems the "twinity" of the signal and idler fields can be shown directly by measuring the 
fluctuations of the difference of the output signal and idler field intensities. 

In our case the signal and idler polaritons do not have the same photon fraction and are 
not coupled in the same way to the external field. Clearly, this should lead to a significant 
reduction of the correlations between the signal and idler output light fields. 



Fluctuations of the output light fields 

Let us first comment on the relevant analysis frequency of the noise. The noise spectra 
vary typically over a range of the order of the polariton linewidth. In noise measurements 
the experimentalists have access to very small analysis frequencies (generally a few tens of 
MHz, i.e. a fraction of fieV) with respect to the polariton linewidths (a few hundreds of 
/ieV). Therefore the noise at zero frequency is the relevant quantity. 

The general expressions of the amplitude noises of the three modes and of the signal-idler 
amplitude correlation can be found in the Appendix. 

In expressions (|57H59|) we have taken equal linewidths for the pump, signal and idler 
polaritons (7*^ = 70 = l2k L )- This assumption is not correct in most microcavity samples. 
Indeed the energy of the polaritons of wave vector 2kt is close to the energy of the nonra- 
diative excitons ; diffusion toward these states is enhanced by their large density of states. 
Moreover, the idler energy is closer to the electron-hole continuum. As a result, the exci- 
tonic linewidth of the idler , ~fb2k L is larger than that of the signal 7^0 and pump 7^ modes. 
The assumption that the cavity linewidth 7^ does not depend on k is correct provided the 
three wave vectors of interest are within the stop-band of the Bragg reflectors. In recent 
experiments, the idler beam has been found to be about 50 times weaker than the signal 
beam (see e.g. Ref. which is consistent with a linewidth ratio 72^/70 = 5. 

We will first give the results in the ideal case (with equal linewidths and an input noise 
equal to the standard quantum noise), and then study the influence of the imbalance between 
signal and idler and the input excitonic noise. 
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Ideal case 



The amplitude noises of the pump, signal and idler beams as well as the signal-idler am- 



plitude correlation are drawn in Fig. 0] as a function of the pump parameter a = yl^/I^ g 
in the case of equal linewidths (and no input excess noise. Although the curves go up to 
cr=5 let us recall that the model is not correct too far above threshold where we can no 
longer neglect multiple diffusions. 

Let us first observe that the signal and idler noise spectra have exactly the same shape. 
However the idler noise is drawn towards the standard quantum level due to its low photon 
fraction which causes important losses at the output of the cavity. It is easy to show that 
the ratio of the noise signals 5 — 1 is equal to the ratio of the photon fractions: 



s * out W-i ci 

A,out,r^\ 1 ~~ Hi \ > 



S*£*(Q) - 1 C. 

The signal and idler amplitude fluctuations diverge close to the threshold (for o — ► 1 + ). 
Noise reduction is obtained above a = 1.55. It grows with the pump intensity and saturates 
at a value .... The amplitudes of the signal and idler beams are very strongly correlated 
slightly above threshold. The correlation tends to one in the vicinity of the threshold (a — > 
1 + ) and vanishes rapidly when increasing the pump intensity. All these results are similar 
to those obtained in nondegenerate OPOs j^g . 



Influence of the signal-idler imbalance 

In this paragraphe we still suppose that there is no input excess noise (n Q = = 
^2k L = 0). Let us compare the results with different linewidths to those of the "balanced" 
case (7 = 7fc L = 7o = l2k L — la — lb) m equations (f70|) - (f7^) . It is easy to show that 
the excess 5 — 1 noises of the pump, signal and idler beams are respectively multiplied by 
la/lk L , la/lo an d l a /l2k L - The signal-idler correlation (without normalization) is multiplied 

by la/ y/l0l2k L - 

As an example the case 70 = ik L = l2k L /5 = la is shown in Fig. El The amplitude 
noises of the pump and signal beams have not been represented since they are unchanged. 
The excess noise and the noise reduction are strongly reduced on the idler beam due to its 
larger losses (Fig. El (a)). The signal-idler correlation remains strong close to threshold but 



17 



decreases more rapidly with increasing pump intensity (Fig. 13(b)). 

Influence of input excess noise 

We have assumed that the biggest source of noise for a given polariton mode is the 
luminescence of an exciton reservoir which is populated by the polariton mode itself. The 
input noise for a given mode is then proportional to the mean exciton number in this mode. 
The efficacy of this process is given by the j3 coefficient introduced above ; here we will assume 
that (3 has the same value for the three modes. Slightly above the oscillation threshold, the 
pump mode is much more populated than the signal and idler population ; then the input 
noise is much greater for the pump than for the signal and idler. 

Fig.lHlshows an example in the "balanced" case for a noise parameter /?=5.10~ 5 , evaluated 
from noise measurements on the light reflected by a microcavity sample |9| . The input excess 
noise cuts down the noise reduction. Its influence increases with the pump intensity since it 
is proportional to the mean exciton population. However the correlation is actually enhanced 
by the excess noise. It is due to the fact that the pump input noise is distributed equally 
between signal and idler and contributes to the correlations. 

The quantum domain 

Our model predicts strong correlations between the signal and idler light fields. When 
can we say that these beams are quantum correlated ? We will use two different criteria, 
one of "quantum twinity" and one associated with QND measurement. 

Quantum twinity 

In degenerate or quasi degenerate OPOs, the signal and idler output beams have the 
same mean field values and the same noise properties. Quantum correlations between them 
are evidenced by measuring the noise of the difference between signal and idler intensities 
and comparing it to the standard quantum level. The idea behind this is to compare the 
fields under consideration to a classical production of twin beams, which can be achieved by 
using a 50% beamsplitter. 
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In our case, one beam is much more intense than the other one (the ratio of the intensities 
is of the order of 10 for equal signal and idler linewidths). What happens if the two light 
fields A% and A 2 under consideration have different mean values and different noises Si 
and S 2 7 To produce classically twin beams of unequal intensities, one can use an unequal 
beamsplitter. The field fluctuations at the output of such a beamsplitter write 



5A X = tSA in + r5A v (63) 
5A 2 = r6A m - t5A v (64) 

with t 7^ r, where A{ n is the input field and SA V the vacuum fluctuations entering through 
the other port of the beamsplitter. Now the difference S A_ = 5 Ax — SA 2 is not helpful, since 
it does not give a quantity which is independent of 5A in . However, the correlation between 
the two beams is independent of SA in : 

(WA 2 ) 2 dassicaltwins = «<L4 2 > - 1) ((SAI) - 1) (65) 
which can also be written as: 

(C c l ass i ca ltwins) 2 = ^1 — ~S~^) \ ~ ~S^) 

We can evaluate the "twinity" of the beams by using the quantity 

G = , l ~ C (67) 

1-J-Ul-i 



Si J V 5 2 

which is a generalization of the usual squeezing factor on the intensity difference. G 
smaller than 1 means that one has been able to produce two fields which are more identical 
than the copies from a beamsplitter. Moreover, it is possible to show that this criterium 
does not depend on the way by which the two classical twins are produced 

Experimentally, one can measure separately C, Si and S2 and compute G from (|67p. One 
can also amplify in a different way the two photocurrents in order to measure the quantity 
5A a = aSAx - 5A 2 /a. When a 2 = y/S 2 /S 1 then: 

G = ^- (68) 
2 VSA- y / (S 1 -l)(S 2 -l) 



19 



G is proportional to the photocurrent fluctuations when the gains are adjusted so that 
the noise levels are identical in the two channels. The denominator in (jUHj) can be evaluated 
from the excess noises of each field. 



QND correlation 

A further level of correlation is achieved when the information extracted from the mea- 
surement of one field provides a QND measurement of the other, so that it is possible, using 
the information on one field, to correct the other from a part of its quantum fluctuations 
and transform it into a squeezed state. This criterium is widely used in the field of QND 



measurement 



. It can be expressed in terms of the conditional variance 



V ll2 = S 1 (l-C 2 ) (69) 

Note that when the two beams have different noises (Si 7^ 5*2) one has two conditional 
variances and therefore two possible criteria. This shows that the QND criterium evaluates 
the correlation from the point of view of one beam, and is not an evaluation of the quantum 
correlation between the two fields. One possibility is to state that the two fields are QND 



Discussion 

We first investigate the " QND criterium" . The conditional variances are shown in Fig. 
in the case of equal linewidths and zero input excess noise. From the point of view of the 
idler beam, the conditional variance is always lower than 1, if only by a few percent. From 
the point of view of the signal beam, the quantum domain is very small: it begins at a = 
1.53, very close to the point where it begins to be squeezed. It is only between a— 1.53 and 
a=1.55 that we get "QND correlations" between beams that individually have excess noise. 
For cr > 1.55 the QND correlation criterium is satisfied, although the correlation is quite 
small, because both beams are squeezed. In conclusion, no significant "QND correlations" 
can be observed on the signal and idler output beams. 

We now investigate the behavior of the quantity G evaluating the "twinity" of the signal 
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and idler beams. It is drawn in Fig. |H] as a function of the pump parameter in various cases. 
In the case of equal linewidths and zero input excess noise, G goes down to 0.85 which 
indicates the "quantum twin" character of the two beams. If we take the nonradiative losses 
of the idler polaritons into account (we set again r y 2 k L — 57o) G only goes under 1 by a 7 
percent. However the input excess noise (corresponding to the resonant luminescence of the 
three polariton modes) has little effect on the quantum correlations. As explained above 
this comes from the fact that the pump input noise (which is the strongest slightly above 
threshold, when the pump polariton population is much larger than the signal and idler 
populations) is equally distributed between the signal and idler modes and helps building 
up correlations. 

In conclusion, in present-day microcavity samples the "quantum twinity" criterium is 
overcome by only a few percent. This is due to the fact that only the polariton fields 
are perfectly correlated, and we can only observe their photonic parts. A simple image is 
the following: we observe the polariton system through a beamsplitter which amplitude 
transmission coefficient is equal to the Hopfield coefficient Co, which leads to losses that 
destroy the quantum effects. The correlations are further reduced by the imbalance between 
signal and idler. The photonic part of the idler is very small (of the order of 0.05) which 
corresponds to large losses. 

CONCLUSION 

We have presented a novel quantum model allowing to calculate the quantum fluctua- 
tions of the beams produced by a semiconductor microcavity in the regime of parametric 
oscillation. It extends the model developed by C. Ciuti et al. above threshold and includes 
the noise coming from the exciton part of the polaritons. 

We show that some quantum correlation exists between the signal and idler beams in the 
vicinity of threshold. Taking the parameters of microcavity samples which have been shown 
to work in the parametric oscillation regime, it can be seen that the correlation overcomes 
the quantum limit by a few percent. The measurement of these correlations would be of 
great interest, since quantum correlations between the output beams, however small, are an 
indication of much bigger correlations between the intracavity polariton fields. For example, 
in the ideal case at threshold (see Fig. ??), if we measure a gemellity G=0.91 this corresponds 
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to perfect correlations inside the cavity. 

In order to observe better quantum correlations between the output beams, it is very- 
important that the signal and idler linewidths should be made as equal as possible. A 
simple solution would be to use a low finesse cavity. Then the nonradiative losses would 
be less important with respect to the radiative losses and the ratio of the signal and idler 
linewidths would be smaller. A compromise has to be found because the oscillation threshold 
would also be higher. 

We acknowledge fruitful discussions with C. Fabre, C. Ciuti, P. Schwendimann and A. 
Quattropani. 



APPENDIX : NOISE AND SIGNAL-IDLER CORRELATION 

In this paragraph we give the general expressions for the amplitude noises of the sig- 
nal, pump and idler output light fields at zero frequency (denoted by S^' out , S^: out and 
Sp 2 ° ut respectively), and the signal-idler amplitude correlation at zero frequency (denoted 

«0 01 



s Aout = 1+C 2 To. 1 



X 



1 + 



A 2 no7 60 7 2fc£ + A 2 2 fcL n 2kL 7fc2fc L 7o 
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1 x [(a - 2) (3a - 2) (X%n olbolkLl2kL (73) 



lk L l0l2k L 

+Xi kL n 2 k^b2k L lk L lo) ~ 8 (cr - 1) X^r2 kL 7 bfcL 7o72fc i ]] 



where no, and n2k L are the input excitonic noises. From these expressions, it is 

A, out 

L 



easy to calculate the normalized signal-idler correlation at zero frequency C a ' o ^ (f2), using 
definition (jo3j) . 
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FIG. 1: Energy dispersion of the two polariton branches for a microcavity sample having a Rabi 
splitting of 2.8 meV, at zero cavity-exciton detuning. The arrows show the parametric conversion 
of the pump polaritons (~ 10°) into signal (0°) and idler (~ 20°) polaritons. 
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FIG. 3: Ratio of the photonic fractions of the signal and idler polaritons as a function of the 
cavity-exciton detuning 5. The Rabi splitting is 2.8 meV. 
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FIG. 4: (a) through (c) : amplitude noises at zero frequency of the pump, signal and idler beams 
respectively, (d) : signal-idler amplitude correlation at zero frequency. The three modes are 
assumed to have the same linewidths, and the input noise is set as equal to the standard quantum 
noise. 
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FIG. 5: (a) Amplitude noise of the idler beam and (b) signal-idler correlation at zero frequency 
as a function of pump intensity for ^2k L = 57o- On both plots, the curve in dashed line is the 
"balanced" case j2k L = 7o- 
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FIG. 6: Noises at zero frequency in the case of equal linewidths, with an input excess noise given 
by [3=(5 c /2. The ideal case /3=0 is represented on each curve as a dashed line, (a) pump beam 
amplitude noise ; (b) signal beam amplitude noise ;(c) idler beam amplitude noise ; (d) signal-idler 
amplitude correlation. 
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FIG. 7: (Dash-dotted line : conditional variance of the signal intensity fluctuations, knowing those 
of the idler ; solid line : conditional variance of the idler intensity fluctuations, knowing those of 
the signal. The dashed line is the standard quantum level 
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FIG. 8: Value of the "gemellity" G as a function of the pump parameter, in three different cases, 
(a) solid line : ideal case where all linewidths are equal and there is no excess noise, (b) dashed 
line : different linewidths for the signal and idler modes j2k L = 5 70, and no excess noise, (c) 
dashed-dotted line : all linewidths equal, and some excess noise given by f3=(5 c /2. 
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